Quantum Spin Current Induced Through Optical Fields 



O 
O 

(N 

Or 



Xiong-Jun Liu a * , Xin Liu b > c , L. C. Kwek a ' d , C. H. Oh Q t and Mo-Lin Ge b 
a. Department of Physics, National University of Singapore, 2 Science Drive 3, Singapore 117542 
Theoretical Physics Division, Nankai Institute of Mathematics, Nankai University, Tianjin 300071, P.R.China 

c. Department of Physics, Texas A&M University, College Station, Texas 77843-4242, USA 
d. National Institute of Education, Nanyang Technological University, 1 Nanyang Walk, Singapore 639798 

(Dated: February 2, 2008) 

We propose a scheme to generate quantum spin current via optical dipole transition process. By 
coupling a three-level system based on the spin states of charged particles (electrons or holes in 
semiconductor) to the angular momentum states of the radiation, we show that a pure quantum 
spin current can be generated. No spin-orbit interaction is needed in this scheme. We also calculate 
the effect of nonmagnetic impurities on the created spin currents and show that the vertex correction 
of the spin hall conductivity in the ladder approximation is exactly zero. 

PACS numbers: 72.25.Hg, 72.25.Fe, 72.20.My, 73.63.Hs 
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Spintronics 0, 0], the science and technology of ma- 
nipulating the spin of the electron for building integrated 
information processing and storage devices, showed great 
promise and developed rapidly in recent years. In prac- 
tical application, one of the most important goals is to 
create spin currents 0, 0] • For this many interestingand 
basic phenomena, e.g. the spin hall effect 01 l1 m 
the system with spin-orbit coupling, have been discov- 
ered and further studied. On the other hand, spin cur- 
rents can also be generated with the interference of two 
optical fields 0, H, @, or with optical Raman scatter- 
ing effects 01 ■ Most of these methods in the generating 
quantum spin currents rely on the spin-orbit coupling 
such as Rashba coupling (caused by asymmetry of quan- 
tum well) and Dresselhaus term (induced by asymmetry 
of bulk crystal), etc. Moreover, coherent manipulation 
of spin states in semiconductor via optical dipole transi- 
tions has theoretically and experimentally been studied 
by many authors fllL [T3 | . With these ideas in mind, we 
propose here a new scheme to generate spin currents. 

In this letter, we consider a realization of spin currents 
via optical dipole transition process in three-level system 
based on the spin states of the electrons (n-doped) or 
holes (p-doped) in semiconductor. The short-range scat- 
tering by the nonmagnetic impurities is discussed and the 
vertex correction of the spin hall conductivity is shown 
to be zero in the present optical model. 

The three-level A-type system that we envisage for the 
scheme is given in Fig. ^ The system is confined to the 
two-dimensional x — y plane with area L x x L y provided 
by the semiconductor quantum well. The particles in the 
system are subjected to an uniform electric field E in 
the +x direction and a +z-directional uniform magnetic 
field (Bo = Boe z ) which can cause energy splitting and 
spin polarization in the z direction (S z = sq, s+, s_, see 



Fig. 1(a)). The transition |s_) — > \sq) is coupled by a 
cr+ light with the Rabi-frequency f2i = f^ " 1 exp(iki ■ r), 
where Cl[ (r) is the slowly spatially varying amplitude 
and ki is the wave-vector. Another er_ light, character- 
ized by the Rabi frequency tt 2 = fi| 0) (r)e^ +k2 ' r ) with 

0,2 (r) the slowly spatially varying amplitude, couples 
the transition |s + ) — > \sq), where <fi = tan~ 1 (y/x) and k2 
is the wave-vector. I indicates that the <7_ photons are 
assumed to have the orbital angular momentum HI along 
the +z direction. As in previous works 0, El [U EH we 




FIG. 1: (a)A-type system based on spin states, (b) An ex- 
ample for dipole transitions in the A-type system provided by 
semiconductor GaAs quantum well, where the excited state 
1 1/2) is a conduction band state, and the two ground states 
correspond to light-hole state (for | — 1/2)) and heavy- hole 
state (for 1 3/2) ) , see ref. |l3|l and references therein, (c) Illus- 
tration of the two-dimensional system for present model. 
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focus on the generation of spin current without particle- 
particle interaction. Defining the flip operators as <t m „ = 
with energy levels /i, v — so,s+,s_, the Hamilto- 
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nian in present system reads 

H = Ho + Hj (1) 

H = -^-{ihVxy + eA ) 2 + V(x) ~ g s fi B B S z 

Hi = -(M2 2 <t SoS+ + hfli& SoS _) + h.c, 

where V ' xy — e x d x + e y d y , V(x) = —eEx is the elec- 
tric potential, parameters e, m e , ps and g s represent the 
charge, the effective mass of a particle, Bohr magneton 
and Lande factor, respectively. Ao is the vector poten- 
tial of the applied magnetic field, i.e. B e z = V x A . 
Note that spin-orbit interaction is not required in our 
model (we consider a very small or zero Rashba term in 
present case, see e.g. |l4|). To facilitate the discussion, 
Hq is written in r-representation in eq Q. The interac- 
tion Hamiltonian Hi can be diagonalized with the local 
unitary transformation: Hi = U(y)HiW(t) with 
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where S — (k 2 — ki) ■ r + l<j> and the mixing angle 
tai\9(r) = |fi 2 /^i |. Under this unitary transforma- 
tion, the three eigenstates of Hi are easily obtained as 

l*a> = Ei=o,± u, ? ) l*J>« where a = The co- 

efficients of these states are = 0, = cos 9, 

= — sin#, Wq^ = "^7f: u '+ : ' = = '-"72 Sm ^ ano - 
— ±-^cos0, and the corresponding eigenvalues 

read = 0,E^ = ±^\Q 2 \ 2 + l^il 2 - Eigenstate |* ) 
is ty pica lly known as a dark state for a three-level A sys- 
tem |15j . 

The Hamiltonian JQ| can be rewritten in a covariant 
form with the definition of covariant derivative opera- 
tor D = — iV — A, where A = Ao + A e is the sum 
of Abelian gauge potential A and the non-Abelian one 
A e = — iU(r)VW(r), and the latter definition does not 
lead to the curvature term since it is a pure gauge. By a 
straightforward calculation, we find the 3x3 matrix A 
has the following form: 



^ sin 2 9VS 
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which generally is a SU (3) non-Abelian gauge poten- 
tial. The off-diagonal matrix elements of it represent the 
transitions between each two of |^o) an d |^±)- How- 
ever, considering the present system is non-degenerate, 
we can introduce the adiabatic approximation when both 
light fields vary sufficiently slowly in space. For a nu- 
merical estimate, we consider the case | ^ | = \p cx p 

where A is constant and p — \/ x 2 + y 2 the distance from 
z axis 0. The typical values can be taken as : 
fi 2 = in 2 + n 2 ) 1 ' 2 ~ lOps" 1 , I < 10 4 , A ~ lO 4 ™- 1 and 
the radius of the interaction region in the x — y plane 
is R ~ 1.0cm, the transition rate can then be evalu- 
ated as r a {, < 10~ 3 <C 1 for any two different states 
(a 7^ b = ±,0). Thus we neglect the off-diagonal ma- 
trix elements in eq.© and yield a nontrivial adiabatic 
gauge connection A — > A dlag , which is associated with 
a 3 x 3 diagonal matrix for the magnetic field with each 
component (B k = ^e kmn F mn = ^e kmn [D m , D n ]) 

B u = B 33 = (B + B e )e z , B 22 = (B - 2B e )e z (4) 

where B e = MX 2 /e indicates the strength of additional 
effective magnetic fields induced by the optical fields de- 
pends on the angular momentum HI that can have a large 
value from a vortex optical beam |l7j . For our present 
purpose we choose Bq = B e /2 so that B\ x = i?| 3 = 



— B%2 = 3MX 2 /2e. In this case, we obtain an interesting 
result that the z component of B 22 is opposite to that of 
Bu (and of B33) (see fig. 2(a) and (b)), a key point for 
realizing quantum spin current. 
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(a) (b) 

FIG. 2: (color online) (a)(b) Schematic illustration of the 
electric field and the effective magnetic field, where B° = 
B 22 ,B+ = B 11 and B = B 33 . 

Up to this point, the interaction Hamiltonian is di- 
agonal and we can rewrite H in eq.Q i n the following 
effective form: 

H eff = ^-D 2 + V(r), (5) 

where V(r) = (U(r)VU(r)) m9 is a 3 x 3 diagonal matrix 
with its elements Vn = V 22 + f2o, V22 = V(x)—g s pBBaS z 
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and Vix = V22 — f^o- The momentum operator is defined 
by p — —iKD and it has the nontrivial communication 
relations: [n,pj] = iHSij, [n,rj] = 0,\pi,Pj] = -iKFij. 
Generally the dynamics of this situation is determined 
by the SU(3) gauge symmetry. By considering adia- 
batic condition, the transitions between any two of the 
states |^ , ct )(a = 0,±) can be neglected and the symme- 
try SU(3) readily reduces to an Abelian one: SU(3) — > 
U(l) ® U(l) ® U(l). The quantum state of each particle 
can be expanded by the complete eigenbasis of the in- 
teraction Hamiltonian |^) = J2 a =o ± ^a\^a), where the 
coefficients can be determined by the initial condi- 
tions. Note that each eigenstate \^ a ) consists of differ- 
ent spin states \s ) and |s±). However, since are 
eigenstates of the interaction Hamiltonian Hj, and H 
does not lead to the transition between \sq) and \s±), 
the expectation value of the z-axis spin polarization of 
any state |\& Q ) is time- independent. Therefore, for con- 
venience, we can treat the states \^ a ) approximately as 
new "effective spin states" with their z-axis spin polar- 
ization calculated by S™ — (^ a \S z \^f a ). It is easy to 



see that S : 



(0) 



cos' t/s + 
« 2 6s- 
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6 s o + \ (sin 2 9s + + cos 2 0s_ ) . Since the particles in differ- 
ent effective spin states \^f a ) experience different effective 
magnetic fields, they may move in opposite directions de- 
pending on the index a, generating a spin current. 

To facilitate the subsequent calculations, we note that 
the particles with effective spin state \^ a )(a = 0,±) 
interact with the magnetic field B Q where B° = B22 
and B ± = Bn = B33. It is convenient to choose the 
diagonal-elements of the gauge connection A(r) as A a = 
-B" X j SO that Py = k is a good quantum number. 
With the definition a ak = -^-[(x- + £+i%)S] 

where B = \B a \ = 3B e /2 and l B = y/hc/eB, the Hamil- 
tonian for a given k can further be rewritten as 

H e ff = hw(al k a ak + - 



gB S« + H c 



(0) 



where H c = — 

eB 



jj^-mc — -gffcc, g = g s fJ-B and uj = 
Operators a a k satisfy the communication relation 



[ a ak, a '/3k'l = fiafidkk'- The eigenstate of above Hamil- 
tonian can be written as \n, k, a) with its eigenvalue 
e n ,fc,a = (n + 1/2) fio; — gB S" + H c . The analytical re- 
sults allow us to calculate the charge and spin current. 
The spin current operator for a single particle is defined 
by if" = %(SzVy + v y S z ), where v y — [y , H] / ih = p y + xuj 
is the velocity operator in the y direction j(| and the 
corresponding charge current operator reads j c — ev y . 

With the above definitions, for a A e -particle system 
the average current density can be calculated by 
1 



J. 



L x 



^ ] Uc~y)nkf( e n,k,a)- 



(7) 



where (jc' y )nk = Y, a ( n i k -> "Ij'cfyK k > a ) is the current 
carried by one particle in the state \n, k). /(e n ,fc,a) is the 



Fermi distribution function, and N e = J2 n k a f( e n.k,a)- 
When the coefficient $ Q of the initial state takes on the 
simple value |$ | 2 = 1/2, \®+\ 2 = cos 2 7/2 and |$_| 2 = 
sin 2 7/2, a pure spin current is obtained by 



ink 



eE x 

mix 2 



cos 20(s-| 



-)-s ), jc = 0. (8) 



The spin current in above equation is dependent on the 
space position (x,y). This is because the effective spin 
polarization of the state \^ a ) is dependent on space (see 
Fig. 3). Together with the Eq. (0 we further calculate 
the average spin current 
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FIG. 3: (color online) Space spreading of generated spin cur- 
rent versus the unit 8 gj *a 1 where A = 10(l/m) for the blue 
line, A = 100(l/m) for the black line and A = 1000(l/m) for 
the red line, p — \J x 1 + y 2 versus the unit m. 
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eE x n e 
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(4s± 



tan-^A^^^l + AV) 
XL xX /l + X 2 y 2 



s± - s ) ,(9) 



where s± = s + — s_, and n e is the filling of charged 
particles (unit I/m 2 ). It is interesting that generally 
the average spin current is still dependent on the po- 
sition y. For practical application, here we consider 
the case that the XL X 3> 1, i.e. the <r_ field is much 
stronger than cr + one. Then from the Eq. 10 we ob- 
tain the persistent spin current J y z — — (s± + Sq) . 
To provide some numerical evaluations, we set 0, ivi\ 
n e = 1.0 x lO w /m 2 JX 2 = 3.0 x 10 13 /m 2 (XL X » 1) and 
the spin polarization s+ = 3S./2, s_ = —h/2, so = h/2 
(fig. 1(b)). For these values, we obtain the spin hall con- 
ductivity a* x = Jy z /E x w 1.38 x I0 2 e. The above result 
can be understood as follow: from the formula of H c 
one can find the particles in the effective spin state |^ a ) 
have a group velocity v Q = —E x /B"e y in the y direc- 
tion. As in present case B + = B~ = —B°, the group 



velocities u H 



Therefore, a pure dissipa- 



tionless spin current can be generated when the particles 
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have the same probability in state \^q) with the sum of 
|* + ) and It should be pointed out that, for GaAs 

quantum well, this numerical evaluation should be re- 
vised by a constant factor, because another independent 
three-level system (s' + = h/2,s'_ = —3h/2,s' = —H/2, 
noted by A' system) can also interact with the optical 
fields besides the one shown in fig. 1(b) [2(J. If we neglect 
the detunings of the optical transitions and according to 
s' ± + s' = 3H/2 < s± + so, the A' system will slightly 
decrease the present generated spin currents. 

By now all the discussions are based on a clean semi- 
conductor system. Noting that the vertex correction may 
remarkably cancel the spin hall currents in Rashba sys- 
tem when the s-wave impurity scattering is considered 
[Isl Hfll ] , it is much suggestive to discuss such effect in 
the present model. For this we consider the following 
short-range nonmagnetic impurity Hamiltonian: 

fli mj ,=5^u«(r-R J ). (10) 

3 
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The free retarded Green's function G R (e.k) — [e — 
Heff}^ 1 is diagonal in the effective spin space: 

GrI k ) = 1 -—-6nn>6 aa >. (11) 

' ' e - e n ,k,a + "7 

In the Born approximation, the self-energy related 
to free green function in our system can be calcu- 
lated by T, R (e) = n irnp u 2 J Ne{e n ,k, a )G { R } na {e, k)de nt k, a , 
where rii mp is the impurity density and A4(en,fe,a) 
the electron state density. It can be verified that 
£ fl (e) = —inni mp u 2 J\f e (e), thus the total Green's func- 
tion reads G R - nia (e,k) = e _ 6n fc 1 a _ H/2T <W<W- Here 
t = 2imi mp u 2 Ne represents the mean scattering free 
time. With these results we calculate the vertex correc- 
tion of charge current j^% tex in the ladder approximation 
as 



n,a n, a v 

+y/nS nt ni + i)G R -n,a(e, k) ) , (12) 



where the advanced Green's function GU;n,a(e, k) — 
{GR- : n t a(e, k)}* . Notice that the Green's function in the 
present system is diagonal in both Landau level and effec- 
tive spin space: G^ R (e,k) oc 5 nn >5 aa ', whereas integral 
function in the right-hand side of above equation is also 
proportional to 5 n , n '+i or 5 n .n'—i, the vertex correction of 
charge current in ladder approximation is exactly zero. 
When we turn to the ladder diagrams of the spin Hall 
conductivity it® , we have to take into account the ver- 
tex correction in the charge current only |18| . Therefore, 
because of this vanishing vertex correction, the spin cur- 
rent obtained in our optical method reproduces its value 
obtained from a clean system without impurities from 
the outset. 

Finally, we present an experimental realization of the 
initial condition required for the result in eq. J7J) . Firstly, 
we can trap the system in dark state using two light 
fields with equal Rabi-frequencies (i.e. = If^j) and 
switch both lights off simultaneously, hence the state 
l*(*o)) = (l s +) — \ s -))/^- in this step no angular mo- 
mentum is needed for the cr_ light field. Secondly, we 
only turn the ct+ light on again. Meanwhile the popu- 
lation in the state |s+) keeps 1/2 unchanged while the 
state |s_) is coupled to the state \sq) by the a+ light 
and the state reads \^{t\)) = (sin|f2i|(ti — t )\s ) + 



cos|fli|(ti — t )\s-) + |s+)/)V2- Thirdly, we adiabati- 
cally turn on the er_ light field with z-directional angular 
momentum and arrive at the adiabatic evolution of the 
state |¥(f)) = j^jge-<^ + ) - ^e*r*|*_) + |* » 

where 7 = f t y/\Sli\ 2 + \Q2\ 2 dt. Obviously, this is the 
state needed for the initial condition. 

In practice, the spreading of light fields have a bound- 
ary in the x-y plane which may lead to modification in the 
Landau energy levels. However, based on the previous 
results 0, this boundary effect can be safely neglected 
when l\ 2 L x L y 3> 1, which is achieved using a light beam 
with a large angular momentum. On the other hand, we 
should emphasize that the particle-particle interaction 
may lead to a renormalization of the Rabi-frequencies 
£li t 2 of the transitions between states \s ) and \s±) cou- 
pled by the light fields 12(1 l21l . However for simplicity, as 
in previous works 0, ELldll^. we do not need to consider 
such effects in the present model. All these interesting as- 
pects will be further discussed in future publications. In 
conclusion, we have proposed and demonstrated a means 
of generating quantum spin current optical dipole transi- 
tion process. The short-range scattering by the nonmag- 
netic impurities are discussed and the vertex correction 
of the spin hall conductivity is shown to be zero in the 
present optical model. 
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